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Abstract 

The 7-ray emission of blazar jets shows a pronounced variabihty and this feature 
provides hmits to the size and to the speed of the emitting region. We study the 
7-ray variabihty of bright blazars using data from the first 18 months of activity 
of the Large Area Telescope on the Fermi Gamma- Ray Space Telescope. From 
the daily light-curves of the blazars characterized by a remarkable activity, we 
firstly determine the minimum variability time-scale, giving an upper limit for 
the size of the emitting region of the sources, assumed to be spheroidal blobs in 
relativistic motion. These regions must be smaller than ^ 10~^ parsec. Another 
interesting time-scale is the duration of the outbursts. We conclude that they 
cannot correspond to radiation produced by a single blob moving relativistically 
along the jet, but they are either the signature of emission from a standing shock 
extracting energy from a modulated jet, or the superposition of a number of 
fiares occurring on a shorter time-scale. We also derive lower limits on the bulk 
Lorentz factor needed to make the emitting region transparent for gamma-rays 
interacting through photon-photon collisions. 
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1. Introduction 

Blazars are active galactic nuclei (AGN) with a relativistic jet pointing at 
a small angle from our line of sight. This feature allows to observe strong 
relativistic effects in blazars emissions, such as the shortening of time intervals, 
and hence rapid variability. Variability at all frequencies is indeed a defining 
characteristics of this class of sources (see for example Ulrich et al. 1997). The 
study of variability at high energies is particularly interesting, because most of 
the bolometric luminosity of these objects is emitted at 7-rays (MeV-GeV-TeV). 

The EGRET experiment onboard the Compton Gamma-Ray Observatory 
already detected rapid variability with time-scales of few hours. As an example, 
EGRET measured in January 1996 a flux increase of a factor of 2.6 in ^ 8 hours 
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in the emission of the weU-known FSRQ 3C 279 (Wehrle et aL 1998). In June 
1995, EGRET observations of PKS 1622-297 led to the detection of a huge flare, 
with a flux increase of a factor of 3.6 occurred in less than 7.1 hours, with a 
corresponding doubhng time of less than 3.8 hours (Mattox ct al. 1997). The 
early results obtained from the Large Area Telescope (LAT) onboard the Fermi 
Gamma-Ray Space Telescope (Fermi) show rapid variability of few hours (see 
Foschini ct al. 2010, Tavccchio et al. 2010, Abdo ct al. 2010). 

Such a rapid variability is useful to constrain the size R of the emitting region 
of blazar jets. The latter is linked to the minimum variability time-scale (tmin) 
by the relationship R < ct„iinS (1 + z)~^, where S is the rclativistic Dopplcr 
factor {6 = [T{l—/3cos9)]~^, where /3c is the emitting region velocity, 9 the angle 
of its motion with respect to our line of sight and F the bulk Lorentz factor). 
Assuming the emitting region to be a sort of spherical blob that increases its 
size following the jet aperture, we can obtain an estimate of its distance from 
the central engine. Calculating the size of the emitting region allows to evaluate 
at which distance from the central black hole the jet originates. This kind of 
result would give important constraints on jet modeling, because understanding 
how the jet originates near the central black hole is nowadays an interesting 
challenge of blazar paradigm (for a closer examination about blazar detailed 
structure see Begelman, Blandford and Rees, 1984). 

An estimate of the mentioned parameters {R, d) results in general agreement 
with the observed timc-scalcs of few hours range. Some results obtained with the 
High Energy Stereoscopic System (HESS), however, give evidence of very fast 
variability at 7-rays in the TeV BL Lac PKS 2155-304 (Aharonian et al. 2007, 
for multiwavclcnght follow-up Foschini et al. 2007 and Aharonian ct al. 2009). 
During an exceptional outburst in July 2006, this source showed an activity 
characterized by a doubling time-scale of ~ 200 s. Another example of very fast 
variability has been detected by the Major Atmospheric Gamma-ray Imaging 
Cerenkov Telescope (MAGIC) during the remarkable activity of the TeV BL Lac 
Mrk 501 in June and July 2005 (Albert et al. 2007). In this case the variability 
time-scale derived during the most active nights was of ^ 5 min. With the 
current values of the parameters mentioned above, results like those obtained 
for PKS 2155-304 and Mrk 501 are about one order-of-magnitude less than what 
one expects, i.e. tmin ~ Rs/c (Begelman et al. 2008). A solution of this problem 
could consist in increasing the value of the Doppler factor, i.e. the relativistic 
beaming of the jet. 

However, very large values of the bulk Lorentz factor could be possible for BL 
Lacs (Begelman et al., 2008), but are problematic in FSRQs. This is because, 
in FSRQs, the high energy emission is thought to be produced through external 
Compton (EC) mechanism, i.e. inverse Compton on a population of seed photons 
external to the jet. The energy density of these external photons is enhanced, 
in the comoving frame, by F^. A too large value of F makes the 7 — 7^6^ 
process important, and 7-rays would be absorbed. This limit is not present in 
BL Lacs, due to their lack of emission lines. 

The Large Area Telescope (Atwood et al. 2009) represents a good instrument 
to improve studies about blazar variability in the MeV-GeV energy band. The 
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great advantage of LAT operating in scanning mode is to offer an efficient mon- 
itoring of tlie whole sky, witli a complete coverage every three hours, that leads 
to the possibility of drawing 3 hours-binned light curves, i.e. the finest binning 
achievable in scanning mode. This featiire can eventually allow to deepen the 
study of minimum variability. Because of Fermi /LAT useful features, our study 
is focused on the analysis of the 7-ray variability of bright blazars detected by 
this telescope during the period from August 4, 2008 to January 25, 2010. 

In the present work we use a cosmology with Hq = 71 km s~^Mpc~^, = 
0.7 and Qm = 0.3. 

2. Data and analysis 

Data of Fermi /LAT have been downloaded from the publicly accessible web 
site of HEASARC^, together with the software of analysis. We focused on five 
sources characterized by a pronounced variability, selected from a public list of 
monitored sources. According to LAT data policy, a source belongs to this list if 
its flux exceeds a minimum value of 2- 10~^ph cm~^s~^ {E > 100 MeV). In order 
to identify the most variable sources, we flrstly made an overall analysis of the 
whole LAT list. From the data published we obtaincid 1 day-binned light curves 
covering all the time period. The minimum statistical confidence level accepted 
for each time bin is TS = 25, where TS is the test statistic described in Mattox 
et al. (1996). This statistical characterization is provided by HEASARC among 
the data published in the monitored sources list. Therefore, according to this 
selection, we generate light curves taking only time bins with a ~ VTS > 5. 
Because of this statistical selection, each source will have a different number of 
time bins. Therefore, we take into account that the less statistically signiflcant 
light curves might have time intervals without reliable data. 

We selected the most variable sources working on the flux deviation from 
the average flux value. In order to characterize this deviation, we calculated the 
xled with respect to a constant flux model centered on the average flux. The 
minimum value accepted is Xred ~ ^- found five sources out of the sample 
that fulfill this condition: 0235-M64, 3C 273, 3C 279, 1510-089, 3C 454.3 (for 
the specific Xred values, see Table 1). We focuss on this smaller group of sources. 

Data of the most variable sources have been analysed by means of the pub- 
lic software package LAT Science Tools v. 9 . 15 . 2, including the calibration 
files (Instrument Response File, IRF P6_V3_DIFFUSE), the maps of the spatial 
diffuse background and the spectrum of the isotropic background, both cos- 
mic and instrumental. The data sets has been flltered to exclude the events 
outside the good-time intervals. The statistical characterization of data sets is 
based on an unbinned likelihood algorithm, implemented in the gtlike task. 
This algorithm has been used to extract the flux and the photon index of the 
source emission, modeled with a single power-law with the form F{E) oc E~^'y 
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Source 


z 


(ii [cm] 


Fx [^Jy] 


ax 


points 




0235+164 


0.940 


1.88 


10^^ 


0.34 


0.44 


138 


11.61 


3C 273 


0.158 


2.33 


1027 


67.77 


0.56 


213 


16.96 


3C 279 


0.536 


9.51 


1027 


0.98 


0.76 


266 


22.80 


1510-089 


0.360 


5.93 


1027 


0.44 


0.40 


263 


92.36 


3C 454.3 


0.839 


l.()9 


10-" 


2.97 


0.70 


404 


18.41 



Table 1: The five sources out of tfie overall sample that we analysed in details, z, cIl, 
Fx and Qj; arc the parameter used to calculate the 5 factor (data from GhiscUini ct 
al. 2010). Npointa is the number of statistically significant points of the light curves. 
Xred represents the deviation of the light curves with respect to a constant flux model 
centered on the average flux {Npoints — 1 is the number of degrees of freedom involved 
in each calculation). 

(where is the photon index). A gross estimation of the spectral parameters 
is done with the optimizer DRMNFB. The (nitput of this first run is used as an 
input for the second optimizer NEWMINUIT. In the following, we will display only 
statistical errors, since we will not make comparisons with other instruments. 
Therefore, we neglected systematic errors, and there is no need of absolute flux 
calibration. Performing this analysis, we obtained for the five selected sources 
the light curves shown in Figure 1. 

3. Minimum Vctriability time-scale analysis 

The first step of our analysis is the determination of the minimum variability 
time-scale, defined as the shortest doubling or halving time among the significant 
flux variations present in the light curves. To ensure the variation significancy, 
we require a difference of 4cr between the flux values of two adjacents bins, where 
a is the largest error bar involved in the variation. Therefore, we calculated the 
equivalent doubling or halving time for each valid variation, in order to define 
the shortest ones as At^in. Table 2 reports the values found for Af„j„, both 
for the rise and decay phases. On average, they are very similar. 

3C 454.3 shows little variability for most of its light curve, becoming more 
active towards the end of the observation period. In December 2009, 3C 454.3 
made an intense outburst, followed by a period of high activity, that reached its 
climax in April 2010, with a second huge outburst. The first event starts with 
the intense activity that can be noticed at the end of our light curve (Figure 
le). These two outbursts are deeply analysed in many works (Ackermann et al. 
2010, BonnoH et al. 2011, Foschini et al. 2011, Tavecchio et al. 2010), where the 
minimum variability time-scale obtained is about 3 hours. More stringent results 
can be found in literature also for other sources, and the respective time-scales 
are included in Table 2. These results are obtained through a deeper analysis 
of high fluxes phases, on which a finer time binning (up to 3 hours) can be 
performed. Hence, the results can be more stringent. We have included in 
Table 2 this more restrictive time-scales and the corresponding values Doppler 
factors (Smin) and blob sizes {Rmax)- 
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Time [MJD-54682(4Aug08)] 
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Figure 1: Daily 7-ray light curves {E > 100 MeV) of the five most variable sources 
selected for the analysis: (a) 0235+164, (b) 3C 273, (c) 3C 279, (d) 1510-089, (e) 
3C 454.3. 
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The obtained time-scales can be used to derive constraints on some features 
of the emitting region. If we assume this region to be a spherical blob in rela- 
tivistic motion, we can estimate a maximum radius, given by the relation: 

1 + z 

where the Doppler factor S is derived from the transparency condition to 7 -rays 
(Dondi, Ghisellini 1994): 



S > 



Atr, 



where F{iyo) = Fx{vx / vq)""- is calculated for both = 1 GcV and i^q = 10 GeV, 
in order to take into consideration two different standard frequencies of incident 
7 photons. The parameter used in this formula for each source are indicated in 
Table 1. Table 2 shows the specific values obtained for every source. 

We found upper limits for minimum variability time scale to be Atmin ~ 1 
day. These results correspond to emitting regions characterized by maximum 
sizes in a range from R ~ 1- 10~^pc ^ 3- 10-^^cm to i? ~ 9- 10~^pc ~ 2.7- lO^^cm 
and by a minimum Doppler factor in the range 5^2 — 7, i.e. consistent with 
results obtained from detailed modeling performed by Ghisellini et al. (2010). 



4. Outbursts time-scale cuialysis 

The second part of our analysis is centered on the interpretation of the 

outbursts. A well defined outbiirst is a substantial variation in Q\ix, lasting at 
least a few days with a coherent trend in rise and decay. We found this kind 
of variation in all the light curves of the selected sources, even if with different 
shapes and intensities. 

The best example of an outburst can be observed during the exceptional 
activity of 3C 454.3 in December 2009, deeply analysed by Bonnoli et al. (2011). 
This outburst lasts 11 days and starts from an already high state (figure 2). In 
the case of 3C 279, the most coherent case of activity is a couple of long-time 
variations in the first half of the light curve, as can be seen from the complete 
light curve in figure 1. The identication of a good outburst in this case is made 
difficult by the large error bars, related to the somewhat low flux of this source. 
In the central part of the 1510-089 light curve we can see three variations that 
can fit our definition of outburst, even if less defined than the one in 3C 454.3. 
3C 273 presents instead a peculiar profile in its outbursts: every substantial 
variation seems to be composed of at least two initial spikes, followed by a more 
gradual decay. Every spike lasts 3 or 4 days, as can be seen in figure 2. 

A first possibility is to ascribe the outburst to the relativistic straight motion 
of a unique emitting region, assumed to be a spherical blob. From the analysis 
of the blob motion feature, we can derive some properties. 
In a time interval of Ate = ti — to, the blob moves with a velocity /3c and covers 
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Results with vo = 1 GeV 
Source Rise Decay 



^imin ^<T ^min Rmax ^imin ^cr ^min ^^-max 

[days] [pc] [days] [pc] 



0235+164 


< 0.8 


4 


4.3 


1.6-10-3 


< 1.0 


6 


4.1 


1.9 


- 10" 


-3 




6 hr* 




5.5 


6.1 - 10-* 














3C 273 


< 0.9 


5 


4.6 


3.2 - 10-3 


< 0.6 


11 


5.1 


2.1 


- 10" 


-3 




5 hr* 




6.2 


9.6- 10-'' 














3C 279 


< 1.5 


4 


3.4 


2.9- 10-3 


< 0.5 


6 


4.1 


1.2 


- 10- 


-3 


1510-089 


< 1.1 


7 


2.5 


1.7- 10-3 


< 0.9 


6 


2.6 


1.5 


- 10- 


-3 




6 hr* 




3.4 


5.4 - 10--* 














3C 454.3 


< 3.5 


8 


4.6 


7.5 - 10-3 


< 3.0 


6 


4.7 


6.6 


- 10- 


-3 




3 hr* 




8.5 


5.0- 10-" 




















(''^Results with vq = 


10 GeV 












Source 






Rise 








Decay 












no- 




^max 




ria 




Rmax 






[days] 






[pc] 


[days] 






[pc] 




0235+164 


< 0.8 


4 


5.3 


1.9- 10-3 


< 1.0 


6 


5.1 


2.3 


- 10- 


-3 




6 hr* 




6.7 


7.5 ■ lO""* 














3C 273 


< 0.9 


5 


6.0 


4.1 - 10^3 


< 0.6 


11 


6.6 


2.7 


- 10" 


-3 




5 hr* 




7.8 


1.2- 10-3 














3C 279 


< 1.5 


4 


4.6 


3.9- 10-3 


< 0.5 


6 


5.6 


1.6 


- 10- 


-3 


1510-089 


< 1.1 


7 


3.0 


2.1 - 10-3 


< 0.9 


6 


3.1 


1.8 


- 10- 


-3 




6 hr* 




4.1 


6.5- 10-" 














3C 454.3 


< 3.5 


8 


6.2 


1.0- 10-2 


< 3.0 


6 


6.4 


8.9 


- 10- 


-3 




3 hr* 




11.5 


6.7-10-" 















Table 2: Results of the minimum variability time-scale analysis: Atmin is the minimum 
variability time-scale, from which we derived the minimum Doppler factor Smin and the 
maximum size of the emitting region Umax- Smin is obtained from the Dondi and Ghisellini 
relation (1994) for the jet transparency to 7-rays, using i/o = 1 GeV^"' and vq = 10 GeV'^'"' 
as standard frequencies of the incident 7-photons. is the number of sigma of the variation 
corresponding to the minimum variability time-scale. * Results obtained from several deeper 
analysis performed during phases of particularly high fluxes of the sources. 0235+164 has been 
analysed by Foschini et al. (2008), 3C 273 by Foschini et al. (2011), 1510-089 by Tavecchio et 
al. (2010), 3C 454.3 by Ackermann et al. (2010), Bonnoli et al. (2011), Foschini et al. (2011), 
Tavecchio et al. (2010). 
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Figure 2: Details of outbursts analysed in 3C 454.3 (a) and 3C 273 (b) 7-ray light curves. 



a distance R = ffcAte- Two photons emitted at the beginning (to) and at the 
end (ti) of this path would be separated by a distance AR = cAtg — RcosO, 
where 9 is the angle of the blob motion with respect to the photons propagation 
direction, i. c. our line of sight. The two photons are received with a time 
separation Atoutburst given by 

^toutburst = —{l + z) = Ate(l - /3cOS^)(l + z) ^ ^(1 + z) 

where F is the bulk Lorentz factor. For blazars, we can assume 6 ~ sin 9 ^ 1/T 
and cos6' ~ F, hence (1 — j3cos9) « F"^. The distance R covered by the blob 
in Atoutburst is: 

F2 

R = /ScAte = (3 cAt outburst-;—, — 

1 + z 

It is likely that during this travel, the blob expands. To estimate the expansion, 
we can use the results of the previous section as a reliable approximation of the 
starting radius Rbo- We assume that the blob expansion follows the jet aperture 
tjj = 10~^ip-i, therefore the blob radius at the end of its motion is: 

/ Rbo \ 

Rbl = IpRstop = V'( — p + At outburst ^ 



cAt ■ T F^ 

» — — \- tPPCAtoutburst— 

1 + Z I + Z 

where we assumed ~ 1/F and hence 5 ~ F. Therefore we can estimate the 
expansion of the blob as: 

-^61 -, , I ^-t^ At outburst -I r> I /D -r^ / outburst / min \ 

R^r'^ "^^^^^ " ^-^ ^ ( — ro — ) 

The blob should clearly expand considerably and change quite dramatically its 
emission properties, contrary to what we observe. 

Even if we introduce a smaller jet aperture ~ 1° ~ 1.7 • 10~^ V-2, of the 
order of the radio component of the jet measured by VLBI (Jorstad et al. 2005), 
we obtain a significant expansion: 

-^61 n n I o { outburst I min \ 

— — w 3.6 V-2 P Fis ' ' 



RbO r . ^ ^ J 

that corresponds to a variation in the volume by a factor 30. Hence, even such 

a thin aperture angle provokes an alteration large enough to change drastically 
the blob emitting conditions. Therefore, such long variations in flux cannot be 
interpreted with a unique blob displacement. For this reason, a different inter- 
pretation is needed. One possibility is that the variable luminosity comes from 
a standing shock. The emitting plasma flows through the shock at relativistic 
velocities, but the emitting region remains the same. Another hypothesis is sug- 
gested by the outbursts of 3C 273 (see figure 2). Indeed, these variations can 
be interpreted as a superposition of single flares with shorter time-scales. This 
model would imply the consecutive formation of blobs and the superposition of 
their emission. 
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5. Conclusions 

The blazar emitting region can be modeled with a homogeneous spherical 
blob placed in the jet, with a size comparable to the cross sectional dimension 
of the jet. Minimum time-scale variability allows to determine the maximum 
size of the emitting region. For the high variability sources that we analysed, 
the time-scales are found to be Af^in ~ 1 day. Performing a different pointing 
with LAT (Foschini et al. 2010), or in cases of very strong outbursts, minimum 
variability could be estimated with more accuracy. The emitting region derived 
from Atmin have maximum sizes that spread from ~ 1 • 10~^ pc to ~ 9 • 10~^ pc 
and lower limits of the bulk Doppler factor in the range 6 ~ 2 — 7. Our minimum 
variability study is completely model-independent, therefore it can be used as 
a constraint on models deductions. 

Using the outburst time-scale, we tried to describe the relativistic motion of 
the emitting region. In the assumed scenario, if the outbursts were produced 
during a straight motion of the blob along the jet, the emitting region should 
expand considerably and change dramatically its emission properties, contrary 
to the observational evidence. Therefore, a different hypothesis is needed, such 
as replacing the formation of a single blob with a more complex mechanism of 
flux modulation. 
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